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Abstract 

In a pioneer work, John Nagle has shown that an Ising chain with 
competing short and long-range interactions displays second and first- 
order phase transitions separated by a tricritical point. More recently, 
it has been claimed that Nagle's model provides an example of the 
inequivalence between canonical and microcanonical calculations. We 
then revisit Nagle's original solution, as well as the usual formulation 
of the problem in a canonical ensemble, which lead to the same re¬ 
sults. Also, in contrast to recent claims, we show that an alternative 
formulation in the microcanonical ensemble, with the adequate choice 
of the fixed thermodynamic extensive variables, leads to equivalent 
thermodynamic results. 


1 Introduction 

In the beginning of the seventies, John Nagle mia analyzed an Ising chain 
with antiferromagnetic interactions between nearest-neighbor sites, and the 
addition of equivalent-neighbor (mean-held) ferromagnetic interactions be¬ 
tween all pairs of sites. Depending on the strength of the competition, this 
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system was shown to display second and first-order phase transitions sepa¬ 
rated by a “special critical point”, which was later named a tricritical point 
[3]. A few years ago, this problem has been revisited by some authors as 
one on the “paradigmatic examples” of the inequivalence of ensembles, in 
which the very localization of the tricritical point was supposed to depend 
on the ensemble (canonical or microcanonical) that was used to carry out the 
statistical calculations Biia- We have recently disproved similar claims of 
inequivalence of ensembles for a long-range version of a spin-1 Ising model [6]. 
In the present article we give arguments to show the equivalence of solutions 
in Nagle's model. 

The Hamiltonian of Nagle's model in zero external held may be written 


as 



( 1 ) 


where dj = ±1 for i = 1, 2, ..., N, the long-range interactions are ferro¬ 
magnetic, Jlr > 0, and the presence of a tricritical point requires antifer¬ 
romagnetic short-range interactions {Jsr < 0). In his original work, Nagle 
obtained exact thermodynamic solutions by two elegant and complementary 
techniques, which already refer to different thermodynamic representations. 
Later, this model was solved by easier manipulations, in the usual canonical 
ensemble [7] |H] [H] • Due to its instructive features, and to a number of miscon¬ 
ceptions in the literature, we begin by reviewing Nagle's solution. We then 
resort to a Gaussian identity to establish the (same) solutions in the usual 
canonical ensemble. Finally, we use the corresponding Ising chain, with the 
exclusion of the long-range terms, to write an entropy function in the micro- 
canonical ensemble. With the appropriate choice of the extensive variables, 
and properly accounting for the long-range interactions, we show that there 
are no discrepancies between canonical and microcanonical results. 

2 Original solutions of Nagle 

In the more detailed solution of the problem, Nagle considers the Hamiltonian 
of an Ising chain with the exclusion of the long-range interactions {Jlr = 0) 
and in the presence of a held H, 


N 


N 



( 2 ) 
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Given the temperature ksT = 1/13 and the held H, we write the usual form 
of the canonical partition function, 


Zj = Zj{T,H,N) = J2^^P 


Wi} 


N 


N 


/3Jsr ^ cTiai+i + /3Hy^ ai 


i=l 


i=l 


(3) 


which can be analytically obtained by the transfer matrix technique. In the 
thermodynamic limit, the associated free energy per site, as a function of T 
and H, is given by 

g,=g,{T,H)~~]r,Z,, (4) 

from which we obtain the magnetization per site. 
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(5) 


We then use a Legendre transformation to write another free energy, fj = 
fi (T, m), which is expressed as a function of temperature T and magnetiza¬ 
tion m, 

fjiT,m)=gjiT,H)+mH. ( 6 ) 

In this one dimensional system, there are no problems of convexity, and the 
held H can be eliminated by using equation (|5]). 

Nagle then remarks that the energy per site of an additional term, of 
mean-held nature, may be written in terms of the magnetization, so that we 
have ^ 

Taking into account that this term depends only on m, the free energy 
/ = /(T,m), associated with Nagle's model, dehned by the Hamiltonian 
of equation ([T]), is given by 

/ = / (T, m) = fi (T, m) - (8) 

This is the central equation of Nagle's treatment. In analogy with a Landau 
expansion, the free energy / (T, m) may be written as a power series in terms 
of the magnetization. 


/ (T, m) = ao (T) -F 02 (T) + 04 (T) -f (T) d-..., (9) 
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from which we obtain the critical line (02 = 0 ; 04 > 0 ) and the location of the 
tricritical point (02 = 04 = 0; Og > 0). It should be noted that this Landau 
expansion is written in terms of a density (the order parameter m) and that 
the coefficients of this expansion depend on the thermodynamic fields (in this 
case, the parameters (3Jsr and (3Jlr)- 

In the Appendix of his article, Nagle mentions an alternative calculation, 
in which the canonical partition function, in zero field, is written as 


N 


Z{T,H = 0,N) = exp 


M=-N 


PJlR 




21 R 


2N 


(fl, S) exp l-fiJsR (iV - 4S)], 

( 10 ) 


s=o 


where R = min {(A^ ± M) / 2 }, S is the number of (+, —) pairs, and 


N / R-l\ / N-R-1 


fN(R^s) ^ ^ ^ 1 ( ^ ^ 


( 11 ) 


Although referring to a future publication, Nagle and Yeo never published 
the combinatorial derivation of {R, S), which corresponds to the number 
of microstates of the system with fixed values of N, M and S (in other words, 
with fixed magnetization m and internal energy u associated with the short- 
range terms). This expression of /at (i?, S) is directly related to the entropy 
in the microcanonical ensemble in terms of the appropriate densities. In the 
thermodynamic limit, the partition function Z (T, if = 0, N) is given by the 
maximum term of the sum over M and S. According to Nagle, all the results 
in zero field have been checked in this alternative formulation, in particular 
the location of the tricritical point. 

The choice of independent variables and the alternative solutions of Nagle 
are already a firm indication of the equivalence of ensembles. The asymptotic 
form of the expression of /at (i?. S'), which is directly related to the entropy 
in the microcanonical ensemble, has been independently written by several 
authors [ 10 ], even in recent work with claims of inequivalence of ensembles 
[1|. The most remarkable deduction has been published by Ernst Ising [11] 
in his famous article of 1925. 
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3 Solution in the canonical ensemble 


The usual canonical partition function associated with Hamiltonian ([T]) is 
given by 


Z = Z {i3Jsr,I3Jlr) = y^exp 

i^i} 

Using the Gaussian identity 


N 


fiJsR E 


i=l 


I^Jlr 

2N 


N 


. 2=1 


( 12 ) 


r'+oo 


we have 


where 


dx exp \—x^ + 2ax] = exp (a^) , 

1/2 /.+ 0 O 




dy exp [-fdNf (?/)], 


f{y) = \jLRy‘^ “ 

and Zj is the canonical partition function of an Ising chain, 

N 


Zi = Zi {^JsR,/dJLRy) = 5]^ exp 

{o-i} 


N 


f^JsR O'iO'i+1 + fdJLRy 


2 = 1 


2 = 1 


(13) 

(14) 

(15) 

(16) 


Also, we remark that these results can be obtained from an application of a 
well-known Bogoliubov identity [8]. 

In the thermodynamic limit we write 


fiy) ~ ^JxRy'^ -^\nX{y), 
where A {y) is the largest eigenvalue of a transfer matrix. 


(17) 


A = exp (/Usr) cosh (/3JiR!/)+[exp (2/?JsR)cosh^ (PJlrv) - 2siiih (2/3Jsr)] 

(18) 

We can analyze the critical behavior from an expansion of the asymptotic 
form of / {y) as a power series in ?/, 


/ (y) — ^0 + ^ 21 /^ + ^^4 y^ + A^y^ + ■■■, 


(19) 
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which is equivalent to Nagle's expansion of the free energy f{T,m), given 
by equation ([2]). The critical line comes from A 2 = 0, with > 0, and the 
tricritical point is located at A 2 = A 4 = 0 , with Aq > 0 . 

If we use Laplace's method to calculate the asymptotic form of the inte¬ 
gral fimi . the saddle-point equation is given by 

_ _ sinh {/SJlrv) [l + cosh {/3JLRy)] 

^ cosh {/3 Jlrv) + 

where 

D = sinh^ WJlrv) + exp {-A/SJsr) , ( 21 ) 

so we have the corresponding free energy per spin, 

<, = 9(T) = -ijiHy-iltiA(50. (22) 

In the next Section we derive again the equation of state fl20l) in the context of 
the microcanonical formulation. As in a typical mean-field calculation, there 
is always a paramagnetic solution, y = 0 , but this solution becomes physically 
unacceptable in the ordered region of the phase diagram. If there are several 
solutions, we have to choose the absolute minima, which corresponds to using 
a Maxwell construction (and to recovering the convexity of the free energy). 
From the equation of state fl2U]) . it is possible to check the location of the 
tricritical point, given by 

I^Jlr = exp {—2/3Jsr) , (23) 


which corresponds to A 2 = 0 , and 

PJhR 


^ ^ exp {2/3 Jsr) - exp {6/3 Jrr) 


= l + exp{2/3JsR), (24) 


which corresponds to A 4 = 0, in agreement with Nagle's findings. 


4 Solution in the microcanonical ensemble 

According to the work of Nagle, it is convenient to begin by considering the 
Hamiltonian of an Ising chain, without the addition of mean-field terms and 
in the presence of an external field H, which can be written as 

N N 

Ui = -JsR Y, -HY^i = U-HM, (25) 

i=l i=l 
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where the energy U refers to the short-range interactions, 


U =-JsR^^CFiC^i+i] M = 


Given U, M, and N, the number of microscopic states associated with 
this system may be formally written as a sum over spin conhgurations of a 
product of two delta functions, 

nj = nj{U,M,N) = (u + JsrJ2 M (27) 

{ai} \ i=l J \ i=l J 

We now introduce integral representations of these delta functions, and use 
the transfer matrix technique to carry out the sums. In the thermodynamic 
limit, it is straightforward to write 


Qi {U, M, N) 


dkidk 2 exp [N f {ki, ^ 2 )], 


where 


/ {ki, ^ 2 ) = kiu - k 2 m -h In A {ki, ^ 2 ), 


A {ki, k 2 ) = exp {kiJsR) {cosh k 2 + [sinh^ k 2 + exp {-AkiJsR)] } , (30) 

with u = U/N and m = M/N. 

The entropy per particle as a function of u and m is given by 


Si = Si {u, m) ~ — In G/ ~ /cb / \^ki, /c 2 j , (31) 

where ki and /c 2 come form the saddle-point equations, {df /dki)j.^ = 0 and 

{df/dk2)^^ = 0 , 

exp 

u + JsR = ^JsR -- — — (32) 

cosh k 2 + 


sinh /c 2 1-1-71 cosh k 2 
cosh k 2 + 
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with 


(34) 


D = sinh^ k 2 + exp {—AkiJsR) ■ 

In the entropy representation, we write the differential form 


dsj 


—du - dm, 

T T 


from which we have the equations of state, 


T~[duJj T~[dmJ 


(35) 


(36) 


It is straightforward to use these equations, together with the saddle point 
equations (l3^ and (l33D . in order to show that 


k2 = I3H-, ki = 13, 


(37) 


which is an evidence of the equivalence of ensembles (in the absence of the 
mean-held interactions). 

We now turn to Nagle's model, with the addition of the long-range terms. 
In the presence of the equivalent-neighbor interactions, the internal energy 
is given by the sum of two terms, 

u = usR + ulr = usR — 2 '^lr (38) 

so that both the energy associated with the short-range interactions, urr^ 
and the magnetization m should be hxed in the microcanonical formulation. 
Therefore, the entropy of Nagle's model is still given by the expression sj, 
as in equation (|3Tll . but with the energy given by equation (138|1 . which leads 
to a new differential form. 


ds 



JLRm 

T 


dm — —dm. 
T 


From this expression we have 


1 

T 



JLRfn _ F 

T T 



(39) 


(40) 


where H^f = H + JiR^n is an effective field, including the external field H 
and the effects of the long-range terms. Thus, in zero external field, we use 
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equations (1^ to write k 2 = /3JLRfn and ki = [3. Inserting these expressions 
into equations (1551) and (1551) . we obtain 

sinh {/3JLRm) [l + ^ 2 ] ,,,, 

™ =- cosh iURnm) + -■ <«> 

where 

D = sinh^ iiSJLRm) + exp {-A/SJsr) , (42) 

which is identical to the equation of state (15(1]) in the canonical ensemble, with 
the identihcation of m with y (and which already leads to the location of the 
tricritical point). In contrast to previous calculations, we do not hnd any 
disagreements in the thermodynamic behavior obtained from calculations in 
different ensembles. 


5 Conclusions 

We revisited the statistical analysis of a spin-1/2 Ising chain with antifer¬ 
romagnetic interactions between nearest-neighbor sites, and the addition of 
equivalent-neighbor ferromagnetic interactions between all pairs of sites. This 
system, which is known to display second and hrst-order phase transitions 
separated by a tricritical point, has been used as one of the paradigmatic 
examples of inequivalence of canonical and microcanonical formulations. In 
contrast to these claims, we give arguments to show the equivalence of ther¬ 
modynamic solutions in different ensembles. 
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